Spectral Singularities and Zero Energy Bound States 
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Single particle scattering around zero energy is re-analysed in view of recent experiments with 
ultra-cold atoms, nano-structures and nuclei far from the stability valley. For non-zero orbital 
angular momentum the low energy scattering cross section exhibits dramatic changes depending 
on the occurrence of either a near resonance or a bound state or the situation in between, that is 
a bound state at zero energy. Such state is singular in that it has an infinite scattering length, 
behaves for the eigenvalues but not for the eigenfunctions as an exceptional point and has no pole in 
the scattering function. These results should be observable whenever the interaction or scattering 
length can be controlled. 
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I. INTRODUCTION 

Experimental techniques are nowadays capable to dis- 
cern subtle phenomena in a great variety of fields of 
physics. Considerable progress has been made in con- 
trolling the interaction of trapped fermions providing a 
unique opportunity for the study of different states of 
the same atomic system under variation of the atomic 
interaction. This is achieved by means of Feshbach res- 
onance techniques [l| using an external magnetic field; 
the method allows changing the scattering length in a 
wide range between negative and positive values. In close 
vicinity of the resonance or a molecular state (a halo 
dimer) the scattering length a is very large and changes 
sign when going from the resonance to the molecular 
state. 

A Bose-Einstein condensate of neutral atoms with in- 
duced electromagnetic attractive (1/r) interaction has 
been discussed recently as another system allowing a tun- 
able interaction [2]. The Gross-Pitaevskii equations de- 
scribing this system at a large critical negative scatter- 
ing length yield the effective absorbing potential. This 
critical value where the onset of the collapse of the con- 
densate occurs could be interpreted as a transition point 
from separate atoms to the formation of molecules or 
clusters 0. In optics, using media with complex refrac- 
tive index [1, @, an abrupt phase transition has been 
demonstrated. It is associated with the appearance of 
non-orthogonal super-modes near spectral singularities, 
so-called exceptional points (EP) [6, 7]. Nano-structure 
devices provide another example, where the presence of a 
quasi-bound state - resonantly interacting with the con- 
tinuum of scattering states - leads to a Fano-Feshbach 
resonance at specific system parameters Q- AH these 
phenomena can be explained by an effective theory of 



open systems described by a non-Hermitian Hamiltonian, 
where complex eigenstates can have square root branch 
points, i.e. EPs. 

We recall that an EP is a singularity of a non- 
Hermitian Hamiltonian where eigenvalues and eigen- 
states coalesce. The properties of EPs were studied in 
earlier experiments using micro- wave cavities [9|. In 
many cases these singularities produce dramatic devia- 
tions from the traditional resonance behaviour, for in- 
stance a particularly strong dependence on the interac- 
tion strength (or scattering length) (see e.g. [HI and 
papers quoted therein) . Even in nuclear physics where an 
experimental manipulation of, say, the scattering length 
appears impossible, an understanding of properties of 
weakly bound nuclei in terms of these singularities is be- 
ing approached [l^ . Also, using a two-state model it was 
demonstrated that a non-Hermitian Hamiltonin can gen- 
erate the binding of unstable states [1^ ; the model mim- 
ics the properties of halo nuclei which form the bound- 
aries for a nuclear valley of stability. A kind of unification 
of nuclear structure and reactions based on an effective 
non-Hermitian Hamiltonian has been suggested to under- 
stand the transition of the nuclear chart from unbound 
to bound limit and vice versa [l4|. 

The focus of the present paper lies on the singular be- 
haviour of the energy eigenstates when the interaction 
(or scattering length) is varied around a bound state at 
zero energy. For / > the energies behave as if the zero 
energy bound state was an EP; yet the eigenfunctions 
and thus the scattering matrix do not share the singular 
behaviour. Nevertheless, the effect upon the cross section 
is rather dramatic, especially for ^ = 1. There the low 
energy cross section 'snaps' from an ^ behaviour to 
a ~ E'^ = finite constant behaviour when a bound state 
and an antibound state coalesce a,t E = (or likewise 



2 



when a resonant state moves to E = 0). 



II. RESUME OF KNOWN FACTS 

It is well known 15] that the single particle scattering 
problem with a radial potential that admits bound states 
can - for angular momentum larger than zero - have res- 
onances near to zero energy. By increasing the (attrac- 
tive) potential strength this resonance state evolves via a 
zero energy eigenstate to a bound state. Actually, both 
the resonance and the bound state appear as a pole of 
the scattering function in the complex fc-plane (fc ~ ^/E 
with E the energy). In fact, a resonance gives rise to 
two poles in the lower fc-plane that are symmetrically 
situated with respect to the imaginary fc-axis. When in- 
creasing the potential strength the two resonance poles in 
the lower fc-plane move toward fc = where they coalesce 
and then continue moving away at right angle in oppo- 
site directions along the imaginary fc-axis. The poles on 
the positive and negative imaginary fc-axis correspond to 
a bound and (usually denoted as) anti-bound state, re- 
spectively. Note that the wave functions associated with 
the resonances and the anti-bound state increase expo- 
nentially at large distance; they are the Gamow states 
[l6j . Although these facts are well known, the aspect of 
the coalescence a,t k = E = for specific potential pa- 
rameters has - to the best of our knowledge - not been 
sufficiently appreciated. We use the term coalescence as 
the merging of the two eigenvalues does not give rise to 
the usual degeneracy being characterised by two indepen- 
dent eigenf unctions. The type of singularity encountered 
here appears to be an EP for the eigenvalues. It therefore 
produces physically dramatic effects. However, as shown 
below, when looking at the eigenstates and the scatter- 
ing function the behaviour is remarkably different from a 
genuine EP as it occurs generically in finite dimensional 
matrix problems [l7| . 

Without specifying the potential the wave functions 
are not known except for their asymptotic behaviour. For 
potentials vanishing faster than at large distances, 
the asymptotic behaviour of the zero energy bound state 
wave function is for orbital angular momentum I 
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(1) 



General statements are known for the behaviour of the 
scattering length, the phase shift and the cross section at 
low energies. For non-zero energy eigenvalues the scat- 
tering length ai is defined by the expansion 



exp2i(5, = 1 + m,fc2'+i + 0(fc^''+^) 



(2) 



where 5i is the scattering phase shift and | exp 2i5i — 
Ip/fc^ is proportional to the cross section. It is further 
known that a; > for the case of a low energy resonance 
and that a; tends to infinity for a zero energy eigenvalue. 
Similarly, it is known that a/ < when a bound state has 
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FIG. 1: Trajectories in the fc-plane of the two resonances 
(horizontal curves) when e < tends to zero. At e = the 
two eigenvalues coalesce at fc = and move for increasing 
e > in opposite directions along the imaginary fc-axis. 



just emerged. In other words, the scattering length has 
a first order pole when an eigenvalue occurs at E = 0. 

What happens to Eq.(2) in this latter situation? We 
now turn to this question. 



III. ANALYTIC TREATMENT 

The motion of the poles of the scattering function, that 
is of the eigenvalues of the Schrodinger equation using 
outgoing wave boundary conditions, seem to indicate the 
typical signature of an EP at i? = fc = 0: a square root 
singularity in the potential strength. Given the potential 
strength vq that produces a bound state at fc = for 
angular momentum / > 0, changing the potential to vo+e 
yields the new eigenvalues as 



(3) 



with finite radius of convergence. The labels 1 and 2 
refer to the two resonances for e < 0; accordingly, for 
e > 0, -El = fc^ and E2 = fc| are the bound and anti- 
bound states energies, respectively. Note that Ei^2 < 
for e > while Ei^2 are complex for e < 0. The occur- 
rence of the square root in (3) is clearly reminiscent of an 
EP and also clearly signals the sprouting out of the two 
energies in different directions depending on the sign of 
e. We stress that the square root behaviour in (3) refers 
to the eigenvalues as a function of the potential strength. 
Associating for the scattering function - as usually - the 
upper/lower fc-plane with the first/second sheet of the 
energy plane, respectively, we clearly see that an emerg- 
ing bound state - a pole of the scattering function in the 
first energy sheet at i?i < - will always have an anti- 
bound state - a pole in the second sheet at i?2 < - as a 
partner. 

We illustrate this result explicitly for a square well of 
width TT and 1 = 1. In this case the p-wave bound state 
at E = k = appears for vq = 1. One finds the explicit 
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FIG. 2: Cross sections versus energy (in arbitrary units) for 
e = (straight line), e < (sharp rising curve) and e > 
(lower curve). 



expansions 



k, = +^(^ + ^, + 0(,3/2)) 
^-(^-i. + 0(.3/^)) 



(4) 
(5) 



indicating nicely the motion of the levels (the poles of 
exp2«(5i) when e changes from small negative to positive 
values as illustrated in Fig.l. 

In a similar way one finds for e = the expansion of 
the scattering function around fc = 



exp2i(5i = 1 - z-vrfc + O(fc^) 
3 



(6) 



while (2) remains valid for e^Q. We note that for Z > 1 
the leading power of e in (4) and (5) persists while that 
of (6) is to be replaced by fc^'^^. 



IV. SINGULARITY AND OBSERVABLE 
CONSEQUENCES 

This last result invokes dramatic consequence for low 
energy scattering. It is illustrated in Fig. 2, where the 
cross section - being proprtional to |1 — exp(2i(5ip/fc^ - 
is drawn for e = (straight line), for e — — 10~^ (sharp 
rising curve) and e = +10~^ (lower curve). Only when 
e = is the low energy cross section a constant in the 
energy (we simply consider the cross section integrated 
over the angles) . As soon as e the cross section starts 
at zero and rises with the quadratic power irrespective 
of the sign of e. This dramatic change at low energy 
scattering should be observable whenever the potential or 
scattering length can be controlled as in the experiments 
discussed in the introduction. The cross sections for the 
different signs of e ^ are almost identical for small 
values of energy. The onset of the sharp rise for e < 
with increasing energy is due to the resonance generated 
by the slightly less attractive potential. This rise depends 
directly upon the magnitude of the negative value of e: 
the smaller |e| the smaller the resonance energy, that is 



the nearer to zero the occurrence of the sharp rise; yet 
the initial rise still remains ^ E"^ . 

This sudden switch in the low energy behaviour is a 
manifestation of the specific singularity for the zero en- 
ergy bound state. Yet, while the spectrum appears to 
have the signature of an EP, neither the wave function 
and hence the scattering function displays the typical 
characteristics of an EP. In fact, while the resonances 
and the bound and anti-bound states are represented by 
a first order pole of exp 2i5i , the two residues (the spectro- 
scopic factors) conspire in the limit e — ?> such that the 
pole does not occur for the bound state at zero energy. 
In other words, the spectroscopic factors vanish when 
the zero energy bound state is approached. This is in 
sharp contrast to the behaviour of a genuine EP in finite 
matrix models or for the coalescence of two resonances 
at a complex energy. There the scattering function (or 
Green's function) has a pole of second order [lol . [l8| . 
The double pole is related to the vanishing norm of the 
coalescing eigenfunctions at the EP (sometimes referred 
to as self-orthogonality fiol]). In other words, there the 
spectroscopic factors become large in close vicinity of an 
EP In the present case, the zero energy bound state 
wave function is normalisable. Therefore, even though 
the energies appear to behave like the levels at an EP 
including the singular behaviour, the eigenfunctions do 
not. Note that an EP cannot appear for a self-adjoint 
Hamiltonian, and the zero energy bound state problem 
falls into this class. 

We emphasise that our findings are expected to be 
valid irrespective of the particular form of the radial po- 
tential as long as it falls off faster than and is less 
singular than at r = 0. In particular, the results are 
valid for an effective single particle problem describing 
nucleon-nucleus scattering. It is here where our results 
are expected to have a bearing even in nuclear physics. 
Note that the bound state at zero energy is "loosely" 
bound as seen by the mild decrease of the wave function. 
As stated above the state has all characteristics of a Fes- 
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FIG. 3: Normalised zero energy bound state wave functions 
for 1 = 1 (dashed) and I = 9 (solid) for potentials of differ- 
ent depths; here the width of the potential has been chosen 
unity implying about a factor three for the quotient of the 
two potential depths. The inset illustrates schematically the 
effective potential. 
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hbach resonance; in fact the slightest change of external 
parameters (scattering length) may turn it into a proper 
resonance or a weakly bound state. Of course, in nuclear 
physics the only change of parameters is moving along 
the isotopic or isotonic line. The situation discussed re- 
sembles that of nuclei on the drip line. In connection with 
halo nuclei we stress that our findings are valid for any 
angular momentum larger than zero. We recall in this 
context that the wave function for a zero energy bound 
state is - for higher partial waves - sharply concentrated 
at the surface of the binding potential, the more so the 
larger the angular momentum. In Fig. 3 we illustrate the 
zero energy bound states for I — 1 and 1 = 9. The inset 
shows schematically the potential indicating the role of 
the centrifugal part of the potential. The thick grey line 
indicates the energy region in which we concentrate our 
discussion of such Feshbach resonance. There the dra- 
matic features for low energy scattering under variation 
of, for instance, the scattering length or the potential 
depth do occur, as has been discussed above in more de- 
tail. 

In passing we note a mathematical subtlety. Even 
though the wave function of the anti-bound state cannot 
be normalised (it grows exponentially like the resonance 
states) , we may calculate the scalar product of the bound 
state and the anti-bound state wave function. The expo- 
nential decay of the bound state wave function outweighs 
the exponential growth of the anti-bound state. As a re- 
sult the integral converges and is exactly zero: the bound 
state is "orthogonal" upon the anti-bound state. It looks 
like a contradiction as the anti-bound- and bound state 
become the more alike the smaller e > 0. The apparent 
contradiction is due to the non-uniform behaviour related 
to the singular behaviour: the exponential growth of the 
anti-bound state is shifted to infinity when e — 0; if the 
limit is taken before integration the scalar product does 
not vanish; it is the integral over the far outside tail that 
brings about the cancellation with the first finite part of 
the integral. Explicit expressions are given in the Ap- 
pendix. 

We suggest that, similar to the Gamow resonance 
states, the anti-bound wave function can have a physi- 
cal meaning. Matrix elements describing transitions con- 
taining resonance states are widely used. In a similar 
vein the anti-bound state may be viewed as a resonance 
at zero frequency with a finite width. With experimental 
techniques available nowadays this could be within reach. 
We recall that the significance of an s-wave anti-bound 
state is known since long for the neutron-neutron system 

M. 



V. CONCLUSION 

To summarise, we have shown that at low energy sin- 
gle particle scattering for angular momentum larger than 
zero there is a dramatic difference for the cross section 
between the special situation of a zero energy bound state 



and the existence of either a resonance or a bound state 
at finite energy. For I = 1 the cross section is constant 
in the case of a zero energy bound state, while it obeys a 
quadratic energy dependence for low energy eigenstates 
with nonzero energy (resonance or bound state). While 
these results have been known in principle pTj . the con- 
nection to a specific singularity has - to the best of our 
knowledge - not been made. In fact, it is remarkable 
that the scattering function has no pole at a resonance 
at zero frequency whereas the cross section remains finite. 
We note that for Z > 1 the pattern for the cross section 
translates into a ^ £;2'-2 behaviour for the zero energy 
bound state and into a ^ i?^' behaviour for the non-zero 
eigenvalues. We think that our findings are relevant, for 
instance, for experiments with trapped ultra-cold atoms 
in either fermion and boson gases, where the scattering 
length can be controlled near the threshold energy. We 
believe that the insights gained also have a bearing for 
the understanding of nuclei around the drip line where 
weakly bound states and larger angular momenta are rel- 
evant. 
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Appendix 

For / = 1 the asymptotic behaviour of the bound state 
wave function is for large distances r 



e-*^''(l + hr) 



and for the antibound state 



1^: 



a— bound/ 



(l-fcar) 



with kb > ka > 0- Obviously the two functions become 
identical for kb — ?► 0. We consider the integral giving the 
scalar product. The tail end section (being the range 
from the zero of |V'a-bound) at r = 1/ka to infinity) reads 



^"'■(l + V) e+'=»'~(l - fc^r) 



r^dr = 



exp(l - ^)kl 



l/ka 



It remains finite in the limit kb ^ since k'^/(kb — ka) 
remains finite (implying ka — 0) as seen from (3); an 
example is given in (4) and (5) when the limits e — are 
taken for kb = —iki and ka = ik2- Note that the integral 
above would vanish if the limit is taken before integra- 
tion. For fca ^ it is the finite value of the tail end inte- 
gral that cancels against the first part of the scalar prod- 
uct integration thus causing "orthogonality" of |'0bound) 
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and I'^a-bound) for ka ^ 0. Yet, when kb = ka = the two functions are identical and have a finite norm. 
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